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We improve the calculations of the -q — > 7r°77 decay within the context of meson chiral 
lagrangians. We use a chiral unitary approach for the meson-meson interaction, thus gener- 
ating the ao (980) resonance and fixing the longstanding sign ambiguity on its contribution. 
This also allows us to calculate the loops with one vector meson exchange, thus removing 
a former source of uncertainty. In addition we ensure the consistency of the approach with 
other processes. First, by using vector meson dominance couplings normalized to agree with 
radiative vector meson decays. And, second, by checking the consistency of the calculations 
with the related 77 — > ir°ri reaction. We find an r\ — > 7r°77 decay width of 0.47 ± 0.10 eV, 
in clear disagreement with published data but in remarkable agreement with the most recent 
measurement. 

PACS: 02.50.+S, 05.60.+W, 72.15.-v 

1 Introduction 

The r] — > 7r°77 decay has attracted much theoretical attention, since Chiral Perturbation Theory 
(ChPT) calculations have sizable uncertainties and produce systematically rates about a factor 
of two smaller than experiment. Within ChPT, the problem stems from the fact that the tree 
level amplitudes, both at 0(p 2 ) and 0(p i ), vanish. The first non-vanishing contribution comes 
at 0(p 4 ), but either from loops involving kaons, largely suppressed due to the kaon masses, 
or from pion loops, again suppressed since they violate G parity and are thus proportional to 
m u — rrid- The first sizable contribution comes at 0(p 6 ) but the coefficients involved are not 
precisely determined. The use of tree level VMD to obtain the 0(p 6 ) chiral coefficients by 
expanding the vector meson propagators, leads to results about a factor of two smaller than the 
"all order" VMD term, which means keeping the full vector meson propagator. All this said it 
has become clear that the strict chiral counting has to be abandoned since the 0(p 6 ) and higher 
orders involved in the full ("all order") VMD results are larger than those of 0(p 4 ). 
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Once the "all order" VMD results is accepted as the dominant mechanism, one cannot forget 
the tree level exchange of other resonances around the 1 GeV region. The ao(980) exchange, 
which was taken into account approximately in [2], was one of the main sources of uncertainty, 
since even the sign of its contribution was unknown. 

After the tree level light resonance exchange has been taken into account, we should consider 
loop diagrams, since meson-meson interaction or rescattering can be rather strong. First of all 
we find the already commented 0(p 4 ) kaon loops from ChPT, but also the meson loops from 
the terms involving the exchange of one resonance. The uncertainty from the latter was roughly 
expected [2] to be about 30% of the full width. 

Another relevant question is that no attempts have been done to check the consistency of 
j] — > 7r°77 results with the related channel 77 — > ir°n. The reason is not surprising since there 
are no hopes within ChPT to reach the a (980) region where there are measurements of the 
77 — > 7r°?7 cross section. On the other hand, the explicit SU(3) breaking already present in the 
radiative vector meson decays has not been taken into account when calculating the VMD tree 
level contributions, and this effect changes the results obtained from VMD estimations by about 
a factor of two. 

The former discussion has set the stage of the problem and the remaining uncertainties that 
allow for further improvement. In recent years, with the advent of unitarization methods, it has 
been possible to extend the results of ChPT to higher energies where the perturbative expansion 
breaks down and to generate resonances up to 1 .2 GeV. In particular these ideas were used to 
describe the 77 — > meson — meson reaction, with good results in all the channels up to energies 
of around 1.2 GeV [3]. With these techniques, and always within the context of meson chiral 
lagrangians, we will address three of the problems stated above: First, the ao (980) contribution, 
second, the evaluation of meson loops from VMD diagrams and, third, the consistency with the 
crossed channel 77 — > n°n. In particular, we will make use of the results in [3], where the 
77 — ► ir°n cross section around the ao(980) resonance was well reproduced using the same 
input as in meson meson scattering, without introducing any extra parameters. 

With these improvements we are then left with a model that includes the "all order" VMD 
and resummed chiral loops. 

2 Mechanisms 

Following [2] we consider the sequential VMD mechanism of Fig. 1 which can be easily derived 
from the VMD Lagrangians involving VVP and Vj couplings 
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where and P are standard SU(3) matrices constructed with the nonet of vector mesons con- 
taining the p, and the nonet of pseudoscalar mesons containing the n, respectively. We also 
assume an ordinary mixing for the <j), the u>, the 77 and 77'. 

From Eq. Q one can obtain the radiative widths for V — > Pj obtaining a fair agreement 
with the experimental data in the PDG but the results can be improved by incorporating SU (3) 
breaking mechanisms. For that purpose, we normalize the couplings so that the branching ra- 
tios agree with experiment. These will be called results with "normalized couplings". In this 
way we are taking into account phenomenologically the corrections to the VP7 vertex from an 
underlying field theory. 

The integrated width obtained using this sequential VMD contribution is T = 0.57 eV (uni- 
versal couplings); T = 0.30 ± 0.06 eV (normalized couplings), where the error has been calcu- 
lated from a Monte Carlo Gaussian sampling of the normalization parameters within the errors 
of the experimental branching ratios. 

Our VMD normalized result is within three standard deviations from the value presently 
given in the PDG: T = 0.84 ± 0.18 eV, but within one sigma of the more recent one presented 
in [4], r = 0.42 ± 0.14eV. There are, however, other contributions that we consider next. 

The contribution of pion loops to ?/ — > 7r°77, evaluated in [2], proceeds, to begin with, 
through the G-parity violating 7/ — > ir°ir + ir~ process but it is proportional to m u — arid, and we 
shall include it in the uncertainties of together with other isospin violating contributions. 

The main meson loop contribution comes from the charged kaon loops at 0(p 4 ) and proceeds 
via 77 — > ir°K + K~ —> 7r°77. Note that these loops are also suppressed due to the large kaon 
masses. That is why the 77 — > 7r°ao(980) — > 7r°77 mechanism was included explicitly, with 
uncertainties in the size and sign of the ao(980) couplings. As commented in the introduction, 
the chiral unitary approach solves this problem by generating dynamically the eto(980) in the 
K + K~ — > 7r°?7 amplitude. 

We can illustrate this approach by revisiting the work done in [3] on the related process 77 — > 
7r°77 where the chiral unitary approach was successfully applied around the ao(980) region. Since 
for the 77 decay the low energy region of 77 — > 7r°?7 is also of interest, we will include next the 
VMD mechanisms also in this reaction. Once we check that we describe correctly 77 — > 7r°77, 
the results can be easily translated to the eta decay. We will finally add other anomalous meson 
loops that are numerically relevant for eta decay but not for 77 — > 7r°?7. 

In [3] it was shown that, within the unitary chiral approach, the 77 — * 7r°?7 amplitude around 
the ao(980) region, diagrammatically represented at one loop in Fig. 2, factorizes as 

— it = {t xK + t AK + K -)t K + K - t7I o v (2) 

with tx+K- ,it°n the full K + K~ — * 7r°77 transition amplitude. 

The first three diagrams correspond to I x k tx+K- ,7r°rj of Eq. Q. The meson meson scatter- 
ing amplitude was evaluated in [6] by summing the Bethe Salpeter (BS) equation with a kernel 
formed from the lowest order meson chiral Lagrangian amplitude and regularizing the loop func- 
tion with a three momentum cut off. Other approaches like the inverse amplitude method or the 
N/D method all give the same results in the meson scalar sector. The BS equation with coupled 
channels can be solved algebraically, leading to the following solution in matrix form 



t(s) = [l-t 2 (s)G(s)}-H 2 ( S ), 



(3) 
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Fig. 2. Diagrams for the chiral loop contribution. 



with s the invariant mass of the two mesons, ti the lowest order chiral amplitude and G(s) a 
diagonal matrix, 6iag(G-^ K , G, ;7r ), accounting for the loop functions of two mesons. These G 
functions were regularized in [6] by means of a cut off. 

In Eq. there is another term, t AK + K -t K + K - y) rr which corresponds to the last two 
diagrams of Fig. 3 where the axial vector meson _K"i (1270) is exchanged. 

In addition to the axial vector meson exchange in loops considered, we have to include the 
loops with vector meson exchange for completeness. In fact, some of the uncertainties estimated 
in [2] were linked to these loops. For consistency, once again we have to sum the series obtained 
by iterating the loops in the four meson vertex. 

Of course, when introducing loops with vector meson exchange we have to consider loops 
involving a K* + or a K*° exchanged between the photons, which were not present at tree level. 

In the -q — > 7r°77 case, the meson loop diagrams correspond to those of Tr r] — > 77 but 
considering the ir° as an outgoing particle. 

Since we are considering all the VMD diagrams and the chiral loops, we still have to take 
into account another kind of loop diagrams [2] which involve two anomalous 7 — * 3M vertices. 
Despite being 0(p 8 ) it has been found [2] that they can have a non negligible effect on the r\ 
decay. 

3 Results 

Using the model described in the previous section, we plot in Fig. 4 the different contributions 
to dT/dMj. We can see that the largest contribution is that of the tree level VMD (long dashed 
line). Let us recall that this is a new result as long as we are using the VMD couplings normalized 
to agree with the vector radiative decays. The resummation of the loops in Fig. 2 using Eq.Q, 
(short dashed line) gives a small contribution (0.011 eV in the total width), but when added 
coherently to the tree level VMD, leads to an increase of 30 % in the 77 decay rate (dashed- 
dotted line). More interestingly, the shape of the 77 invariant mass distribution is appreciably 
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Fig. 3. Contributions to the two photon invariant mass distribution. From bottom to top, short dashed line: 
chiral loops from Eq.j2j; long dashed line: only tree level VMD; dashed-dotted line: coherent sum of the 
previous mechanisms; double dashed-dotted line: idem but adding the resummed VMD loops; continuous 
line: idem but adding the anomalous terms, which is the full model presented in this work (we are also 
showing as a dotted line the full model but substituting the full t K +K - jJ?7r o amplitude by its lowest order). 



changed with respect to the tree level VMD, developing a peak at high invariant masses. The 
resummed VMD loops leads, through interference, to a moderate increase of the r\ decay rate 
(double dashed dotted line), smaller than that of the chiral loops considered before. The last 
ingredient is the contribution of the anomalous mechanisms (continuous line), leading again to 
a moderate increase of the rj decay rate, also smaller than the chiral loops from Eq.(|2j. These 
anomalous mechanisms have a very similar shape to the tree level VMD and interfere with it in 
the whole range of invariant masses. Altogether the final result is 



r(r) -> 7T°77) = 0.47 ±0.10 eV (4) 

where the theoretical error have been obtained considering the uncertainties from the vector 
meson radiative decays, the contribution of the l" 1 axial-vector mesons and the isospin violating 
terms. 

Note that although we have considered a new error source from the uncertainties in the vector 
radiative decays, which turns out to be the largest one, we still have reduced the uncertainty from 
previous calculations. 

The result of Eq. is in remarkable agreement with the latest experimental number T = 
0.42 ± 0.09 ± 0.08 eV [4], based on 1600 events and also with those of [5] of T = 0.35 ± 0.20 
based on a smaller statistics of about 120 events, and lie within two sigmas from the earlier one 
in the PDG T = 0.84 ± 0.18 eV. Confirmation of those preliminary results would therefore be 
important to test the consistency of this new approach. Furthermore, precise measurements of 
the 77 invariant mass distributions would be of much help given the differences found with and 
without loop contributions. 
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Finally, we would like to make some precisions concerning the comparison of recent ex- 
perimental data with former calculations of chiral perturbation theory plus VMD estimates. The 
chiral models prior to our work have to be looked in perspective. They rely upon data on radiative 
decay of vector mesons which has changed considerably in recent years. We use updated data 
and find that if the earlier calculations would have used the present data they would get eta decay 
widths about one half what they got. Thus, comparing our result and the old ones is somewhat 
misleading. And it is also misleading to compare the experiments with these calculations without 
this warning. We find also inappropriate to compare experimental results with what chiral per- 
turbation theory would give at order 0(p 6 ) since, as we have mentioned here, these results are 
obtained from a VDM model projecting over 0(p 6 ) the results of the full model, which provides 
a width about a factor two larger than its 0(p 6 ) projection. In other words, a strict chiral per- 
turbative calculation in terms of powers of momentum is not practical for this problem, and the 
full model predictions, which accounts also for higher powers of momentum, have to be taken 
for reference. 

In addition, previous works in the literature have large uncertainties from ignorance of the 
ao (980) contribution. We have improved all this and other small things and evaluate theoretical 
errors from different sources. In summary the most relevant things of the present work are: 

1) Use of updated data for radiative decay of vector mesons. If previous chiral calculations 
are updated making use of these new data the result comes about a factor of two smaller. 

2) Unitarized chiral perturbation theory allows to determine precisely the contribution of the 
ao (980) resonance, which was a major source of uncertainty in the past. 

3) A careful determination of theoretical errors has been done, which was also lacking before. 
Our work has made the most accurate calculation, so far, withing the framework of chiral 

dynamics and can be taken as reference of what chiral theory predicts for this reaction. Previous 
chiral works have to be seen in their value as making the first predictions for this ratio, but should 
not be used as a reference for a quantitative prediction since they contain intrinsic uncertainties 
of more than a factor two. 
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